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ABSTRACT 

We analyse 4.5, 8 and 24 /im band Spitzer images of six gamma ray burst host galaxies at redshifts 
close to 1. We constrain their star formation rates (SFR) based on the entire available spectral energy 
distribution rather than the 24 /im band only. Further, we estimate their stellar masses (M*) based on 
rest frame K band luminosities. Our sample spans a wide range of galaxy properties: derived SFRs 
range from less than 10 to a few hundred solar masses per year; values of M* range from 10 9 to 10 10 
M Q with a median of 5.6 x 10 9 Mq. Comparing the specific star formation rate (4> = SFR/ill*) of our 
sample as a function of M* to other representative types of galaxies (distant red galaxies, Lya emitters, 
Lyman break galaxies, submillimeter galaxies and z ~ 2 galaxies from the Great Observatories Origins 
Deep Survey-North field), we find that gamma ray burst hosts are among those with the highest qb. 

Subject headings: cosmology: observations — dust, extinction — galaxies: high redshift — galaxies: 
ISM — gamma rays: bursts — infrared: galaxies 
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1. INTRODUCTION 

A canonical model is well established for long dura- 
tion gamma ray bursts (GRB): association with stel- 
lar core collaps e events and hence with high mass star 
formation (e.g. iHiorth et alJ 12001 IStanek et~aTl 120031 : 
IZeh. Klose k Hartmannl 12004 ICampana et all I2006T ). 
The emerging picture, however, is co mplex. Most 
GRB host galaxies are fai nt and blue (|Fruchter et al.l 
119991 iLe Floc'h et~aTl I2003T ) . A few hosts show tenta- 
tive evidence of very high sta r formation rates ( SFRs ; 
ICharv. Becklin fc Armusri2002l : iBerger et al]|2003t) . but 
their optical properties do not appear typical of the 
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Camera: iFazio et al.ll2004ft and MIPS (Mult iband Imager 
Photometer for Svitzer; iRieke et al.l l2004h photometry, 
together with optical, near IR, submillimeter and radio 
observations, can help establish how GRB hosts relate to 
other high redshift galaxy populations. This is essential 
if we are to understand the full range of properties of 
star forming galaxies at high redshifts and exploit the 
potential of GRBs as more general probes of cosmic star 
formation. 

In this Letter we study a su bsample of the 16 GRB 
hosts observed with Spitzer bv ILe Floc'h et all (2006). 
We compute SFRs, dust masses (Md us t) and stellar 
masses (M*) for hosts at redshifts close to 1. This red- 
shift is particularly relevant bec ause it has been ar gued 
that the global SFR peaks there (|Madau et al.lll998f) . To 
determine SFRs we fit full (optical to radio) spectral en- 
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ergy distributions (SED), which allows us to calculate 
the most accurate and robust values for the SFR in a 
sample of GRB host galaxies to date. At z ~ 1 -Mdust is 
constrained by 850 /im SCUBA observations. To deter- 
mine M* we fit observed 4.5 p m fluxes. The M* e stima- 
tor is well calibrated for z ~ 1 (|Labbe et al.l 12005) since 
the 4.5 /im observed wavelength corresponds to the rest 
frame K band. This gives us, for the first time, an accu- 
rate value of the stellar mass in those hosts. We assume 
an Vt m = 0.3, — 0.7 cosmology with Hq = 70 km s" 1 
Mpc" 1 . 

2. DATA 

From the sample of 16 GRB hosts (GTO programme 
76) we selected six that have z ~ 1 (lio^)- The reader 
is referred to ILe Floc'h et alj ([2006) for a detailed de- 
scription of the full data set. Each host has been imaged 
with IRAC and MIPS. IRAC observations were 4.5 /im 
(300s per host; scale = 1.220" pixel" 1 ; FoV = 5.21 x 
5.21arcmin 2 ; instrumental PSF FWHM = 1.98") and 
8.0 jum (300s per host; scale = 1.213" pixel" 1 ; FoV = 
5.18 x 5.18arcmin 2 ; instrumental PSF FWHM = 1.72"). 
MIPS observations were 24 /im (420 s per host; scale = 
2.45" pixel" 1 ; FoV = 5.23 x 5.23arcmin 2 ; instrumental 
PSF FWHM - 6"). 

We used official Spitzer Post Basic Calibrated Data 
products (carefully verified with our own reductions). 
Host extraction was based on archival imagery world co- 
ordinate system calibration and visually confirmed with 
optical comparison images from the literature. The 
median separation between the host centroid in each 
Spitzer image and the best set of coordinates published 
was about 1". We measured the flux densities over a cir- 
cled area of radius 2 pixels in IRAC and 3 pixels in MIPS. 
Aperture corrections were then applied to account for the 
extended size of the PSF. Our photometry , prese nted in 
Table [U is consistent with lLe Floc'h et all (|2006h . 

3. SPECTRAL ENERGY DISTRIBUTION FITS 

Photometry available for our sample includes our 
Spitzer flux densities and data points from the litera- 
ture, spanning from optical to radio wavelengths. We 



2 Castro Ceron et al. 



TABLE 1 

Hosts: flux densities, star formation rates, dust and stellar masses 





Redshift 
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GRB Host 
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(i*Jy) 


(i*Jy) 
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-^8-1000 


(10 7 M Q ) 


(1O 9 M ) 
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0.83 


_1 


<3.1 
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2.5 C 


2 


0.4-26 


0.3-1.1 


1.5 


970828 


0.96 
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3.9 ± 1.1 
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94 ± 17 


l.l d 
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30 ± 8 


1.3 


2.5 


980613 


1.10 


3 


37 ± 1 


33 ± 8 


169 ± 36 


70 c 
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428 ± 51 


19 


31 


980703 


0.97 


4 


11 ± 2 


<24 


<85 


30 c 


7 


3.8-226 


2.7-10 


7.2 


981226 


1.11 


5 


4.5 ± 1.4 


<31 


<87 


1.2 d 


5 


1.0-84 


0.7-3.7 


3.9 


990705 


0.84 


6 


19 ± 1 


<18 


159 ± 31 


~5 d 


6 


4.5-173 


3.2-7.7 


9.2 



Re fer ences. — JlTllBloom et al.l|1998l: (H) IDiorgovski et al.ll2001l: <t3h IDiorgovski. Bloom fc Kulkarni 
20031 : ph IDiorgovski et al.lll998l : (|5l) Ichristensen et al.ll2005l : ijai TLe Floc'h et alj|2002l : (|7t) ICharv et al 
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Upper limits are quoted at the 3a level, while errors are la. h We have corrected all the SFRs 
quoted to ac count for diffe rences in their IMF scales with respect to our choice of a Salpeter IMF 
(0.1-100 M ;[Salpctcr 1955). c Corrected using the /3 slope technique (Charv et al. 2002, and references 
therein), typically large r than the Balmer lines decrement correction. d Uncorrected for internal extinction; 
Christcn sen et al.l l2u05 argue for no extinction in the host of GRB 981226. 



scaled a set of SED templates to fit these data: Arp 220, 
an archetypal ultraluminous IR galaxy (ULIRG); NGC 
6946, a well character ised, blue, star forming galaxy 
(both ISilva et alj I1998Q : and 64 SEP templates rang- 
ing from starbursts to quiescent ones (|Dale et all 120011: 
iDale fe Heloul l2002T >. The spectral width of the differ- 
ent bandpass filters was taken into account. Fitting was 
evaluated by means of a weighted least squares method. 

Figure [T] shows our models. The Arp 220 template best 
fits the hosts of GRBs 970828 and 980613. For these two 
hosts the data rule out the other SED templates; we have 
plotted the NGC 6946 template for comparison. For the 
hosts of GRBs 970508 and 981226 we could not unam- 
biguously discriminate a best fitting SED template, so 
we plotted those models yielding the highest and lowest 
values for the SFR (see §[Oand Table [J). Model degen- 
eracy was expected in these two cases because we only 
have optical-TR data. For the remaining two hosts we 
also plotted those models yielding the highest and lowest 
values for the SFR. For GRB 980703, the Arp 220 tem- 
plate approximately reproduced the radio flux densities 
but was inconsistent with the 24 /im upper limit, while 
the NGC 6946 template was consistent with the 24 fim 
upper limit but underestimated the radio flux densities. 
For GRB 990705, the Arp 220 template overestimated 
the 24 /im flux density and was marginally inconsistent 
with the 8 /im upper limit, while the NGC 6946 tem- 
plate underestimated the 24 /im flux density. Reproduc- 
ing the SEDs of these two galaxies is problematic and 
may require dust w ith propert i es diff erent from those in 
our te mplates. The lDale etaTI (|2001h and lDale fe Heloul 
(2002) templates yielded SFR values between those of 
the Arp 220 and NGC 6946 models for all the hosts in 
our sample. 

For each host we may be fitting subcomponents tha t 
differ in their properties (C harmandaris et al.l l2004f ). 
For instance, for GRB 980613 so me components de- 
tected in the optical/near IR bands (jHjorth et al.ll2002l 
IDiorgovski. Bloom fe Kulkarnill2003D are offset b y >2.5" 
from the Svitzer centroid (|Le Floc'h et al.|[2006T) . While 
such effects might induce some scatter in the SEDs we 
can still utilise the SED templates as powerful diagnostic 
tools for the SFRs of the sample. 



3.1. Star Formation Rates 

Using our SED fitting models, we calculated the SFR 
for each host using IR luminosities. We converted flux 
densities into luminosity densities using L v (v IC& t) — 
±7rDlf u (v obseTved )/(l + z) (|Hogg et al.l l20f)l). where D L 
is the luminosity distance, and scaled the SED templates 
to match the data points. IR luminosities (Lg-iooo) 
were obtained integrating under the scaled SED tem- 
plates from 8 to 1000 /im (rest frame). This wavelength 
range was chosen so the SFRs could be co mputed us- 
ing S FR(M Q yr -1 ) = 4.5 x lQ~^L m ergs" 1 dKennicuttl 
1998). The results are summarised in Table Q] Er- 
rors quoted are statistical and assume the template is 
a good representation of the data. In addition there 
may be significant systematic errors related, for ex- 
ampl e, to how well the template represents the actual 
SEP (lMichalowski| [2006l). the L$ 

-looo to SFR conversion 
(jKennicutfj Il998f ) (both of the same order, ~30%), or 
a fact or of two from th e choice of initial mass function 
(IMF; lErb et al.ll2006al) . 

For the hosts of GRBs 970508, 980703, 981226 and 
990705, the lower end of their SFR ranges indicates low 
star formation, consistent with the estimates from the 
UV contin uum and, in GRB 981226, with no in ternal 
extinction (jChristensen. Hjorth fc Gorosab el 2005) . The 
host of GRB 970828 is a moderately star form ing galaxy, 
in good agreement with lLe Floc'h et al. (2006). The host 
of GRB 980613 is characterised by high star formation 
activity. Our SFR value f or this host is ~ 5 time s higher 
than the one obtained by iLe Floc'h et al.l (|2006h with a 
lower uncertainty, because we have fitted the entire SED, 
as opposed to only the flux density at 24 /im. We verified 
that if we base our calculations exclusiv ely on the 24 u,m 
flux d ensities we reproduce the results in lLe Floc'h et all 
(|2006h for all hosts. 

3.2. Dust Masses 

Dust emission dominates submillimctcr wavelengths. 
The total Afdust in a galaxy can be estimated from its 
rest frame 450 /im flux density: Md ust = S v D\j{\ + 
z)K,(y)B{v 1 T)] where S u is the flux density at an ob- 
served wavelength corresponding to the rest frame wave- 
length of 450 /im at z = 1, interpolated from the fitted 
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Fig. 1. — SED fits of the six GRB host galaxies. Filled squares: 
detections (error bars within the squares). Arrows: 3 cr upper lim- 
its (value marked by the base of the arrow). For each host in our 
sample we plotted the fit to the template that yielded the highest 
(Arp 220) and the lowest (NGC 6949) SFR value. Solid lines: best 
fit. Dashed lines: fits we could not discriminate between. Dotted 
lines: fits inconsistent with the data, sh own for illustrative pur - 
poses. Datap oi nts: o p tical/near IR from IVreeswi ik et al .l 1119991). 
ISokolov et ajj ll20Ull). IDioreoyski et al.l 1 I2OOII) , ILe Floc'h et at] 
(120021 ) and idhristensen et al.l (120051 ). Ther mal IR from our pho - 
tometry of Spitzer a rchival data. Far IR f romlHanlon et al.[ (2000). 
Submillimeter from ISmith et al.l 1 119991) , ITanvir et al.l (120041) and 
our ph o tometry of SCUBA arch i val data. Ra dio from Brem er et al.l 
(fl998h. IShepherd et all (IT998I). iBereer. Kulkarni fc Frail (120011 ). 
IDiorgovski eta l. (2001) and Bcrgcr et al. (200^1 



SED templates; v is the frequency (666.21 GHz) corre- 
sponding to a wavelength of 450 /im; k oc v@ is the mass 
absorption coefficient with (3 being the du st emissivity 
index ; and B(^, T) is the Planck function (| Taylor et al.l 
2005). We assumed optically thin dust emitting a grey 
spectrum. This me thod yields statistical errors of ~25% 
(jTavlor et al.1 12005). The ranges in Md us t are estimates 
of the systematic, model dependent error. Our results 
are listed in Table [T] The derived median, Md us t = 8 x 
10 6 Mq (calculated from the lowest Md us t value for each 
host), is consistent with the distribution for starburst 
galaxi es (4 x 10 5 to 7 x 10 8 M©) found bv lTavlor et all 
(|2005l) . 

4. STELLAR MASSES 

W e estimated M* from rest frame K band fluxes 
(e.g. iGlazebrook et al.ll2004fh M/Lpc depends to some 
extent on the composition of the stellar population 
(|Portinari. Sommer-Larsen fc Tantalol l2004h or, accord- 
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Fig. 2. — if as a function of M* for our GRB host galaxy sample 
and other repres entative types of galaxies (for a similar plot see 
Erb et al. 2006b). Black squares: SFR values constrained with a 
best fit model. Grey squares: highest/lowest SFR values for those 
hosts for which a best fit model could not be established. Dashed 
diagonals: SFRs of 1000, 100 and 10 Mqjt -1 respectively. Right 
axis: SFR timescale (Tgpjt = M*/SFR), the inverse of <fi. On this 
scale the solid horizontals represent the age of the universe for the 
marked redshift. Our six hosts clearly have Tgppt < t un j versc , al- 
lowing for a history of constant star formation. The distribution of 
our sample in parameter space suggests that GRBs trace galaxies 
that are not selected with oth er techniques. Data point s: GRBs 
from this Letter. DRGs from Ivan Dokkum etall pool) . LAEs 
from lGawiser et al.l (120061) . The point plotted represents an aver- 
age value of the LAE population as a whole, obtained from stacked 
photometry , and the dott ed ellipse its uncerta inty. Spectroscopi- 
cally c onfirmed LBGs from Shaplcy et al. (2001) and Barmbv_et_al 



(p004Tl . SMGs from IChapman et al.l 120051 m!) and IBorvs eTal 
(2005, S FR), where we ha ve considered Lbol — -LfarlRi then a] 
plie d thelKennicuttl 111998) calibration, z ~ 2 from lErb etaT] (1200: 
and lReddv et al.l (12006V 



ing to lLabbe et al.1 (|2005f ) who used iBruzual fc Charlotl 
(2003) with a Salpeter IMF, on the rest frame U - 
V colour, age and M*. GRB hos ts are blue, young 
and faint (e.g. ILe Floc'h et al.l l2003l iBerger et aH l2003l: 
IChristensen et al.l I2004D; thus, we assum e M/Lk ~ 
0.1 (lowest detection in lLabbe et al.l [20051 is 0.16), ob- 
taining a robu st lower limit. Table ffl summarises our 
M* estimates. IVan der Wei et all (|2006f ) examined red- 
shift dependent systematics in determining M* from 
broad band SEDs. They found no significant bias for 
IBruzual k. Charlotl (|2003l ) models with a Salpeter IMF 
(|Salpeterlll955l ). 

5. DISCUSSION 

We have found the hosts in our sample to span a wide 
range of properties. Their SEDs are fitted with templates 
that vary from a blue, star forming galaxy to a ULIRG. 
Their SFRs and M* are quite different, ranging from the 
host of GRB 980613, which is forming a few hundred 
solar masses a year with M* = 3 x 10 10 M , to the host 
of GRB 970508, which is forming of the order of 10 Mq 
a year with M* = 1.5 x 10 9 M Q . 

We find that our SFR values are significantly higher 
(up to a factor of 30) than the lower limits from the 
rest frame UV continuum emission (L\jy; see Table [1} 
and higher (by a fact or of 6) than thos e corrected with 
the (3 slope technique ([Chary et al.ll2002D . GRBs 970828, 
980613 and 980703 illustrate how even the best estimates 
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of dust extinction in a galaxy from the UV slope ma y fall 
short, not only for ULIRGs (|Charv k. ElbazlboOlD . but 
also for LIRGs (i.e. 10 12 Lq > L 8 -iooo > 10 11 L e ). Our 
SFR value for GRB 981226 is consi stent with no internal 
extinction (jChristensen et al.ll2005h . 

The specific star formation rate (4> = SFR/M*) gives 
an indication of how intensely star forming a galaxy is. 
In Figure[2]we have plotted 4> versus M* for ours and five 
other representative galaxy samples: distant red galax- 
ies (DRG), Lya emitters (LAE), Lyman break galaxies 
(LBG), submillimitre galaxies (SMG) and an ensemble 
of optically selected, z ~ 2 galaxies from the Great Ob- 
servatories Origins Deep survey-North field. GRB hosts 
have som e of the highest 4> value s, as previously sug- 
gested by iChristensen et al.l (|2004[ ). represents the 
SFR timescale, so high values are indicative of GRBs 
tracing young, starbursting galaxies. 

The different methods that have been used to derive 
the SFRs of the various samples plotted in Figure [2] likely 
introduce systematic offsets. Likewise, to determine M+, 
model fitting to the full SED could better account for 
variations in M/Lk- For high SFR galaxies, ongoing star 
formation contributes to the near IR emission, especially 
as they lack an old stellar population. High <fi galaxies 
are particularly vulnerable to this effect. Nevertheless, 
Figure illustrates where, within the larger picture, our 
hosts fall. 

We have shown the capabilities of IR observations to 
characterise GRB host galaxies and have compared val- 
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